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Abstract 

We calculate the exact kinetic evolution of cosmic neutrinos until complete de- 
coupling, in the case when a large neutrino asymmetry exists. While not excluded 
by present observations, this large asymmetry can have relevant cosmological conse- 
quences and in particular may be helpful in reconciling Primordial Nucleosynthesis 
with a high baryon density as suggested by the most recent observations of the 
Cosmic Microwave Background Radiation. By solving numerically the Boltzmann 
kinetic equations for the neutrino distribution functions, we find the momentum- 
dependent corrections to the equilibrium spectra and briefly discuss their phe- 
nomenological implications. 



PACS numbers: 14.60. St, 26.35,+c, 95.35.+d 
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1 Introduction 

It is a common point of view to treat as instantaneous several (prolonged) processes that 
occurred in the Universe during the first instants after the Big Bang. One the most 
striking examples is provided by reheating after inflation. In this case, the approximation 
of treating all inflaton quanta as decaying at the same time 1/r (with T denoting the 
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inflaton decay rate) is useful for many purposes, such as the computation of the reheat- 
ing temperature T R , but it may fail to reproduce the abundances of the relic particles, 
especially if they have masses larger than Tr and have never been in equilibrium with the 
thermal bath Jl|. 

Even when we consider the more common situation of a particle which was first in 
equilibrium with the radiation and then decouples, a detailed study of the whole pro- 
cess may provide more accurate results with respect to the standard computation, which 
assumes instantaneous decoupling. A particularly interesting case is that of neutrinos, 
which were coupled to the thermal bath via the well known weak interactions. 

Let us first consider the standard situation, where neutrinos are assumed to decouple 
when the temperature of the primordial universe is about 2 — 3 MeV [[| , namely when the 
rates of the weak interactions which couple them to the electromagnetic plasma become 
smaller than the Hubble parameter. After this instantaneous decoupling, the neutrino 
spectra maintain their equilibrium shape, with the temperature redshifted as the inverse 
scale factor, T v ~ 1/R(t), due to the expansion of the Universe. In the meantime, 
also the temperature of the electromagnetic bath scales in the same way, until it reaches 
the electron mass T 7 = m e ~ 0.5 MeV. At this stage the e + e~ —>■ 77 annihilations 
occur, without affecting the relic neutrinos previously decoupled. As a consequence, 
the temperature of photons increases with respect to T v until it reaches the well-known 
asymptotic ratio T 7 /T u = (11/4) 1/3 ~ 1.401. 

Relaxing the assumption of instantaneous neutrino decoupling, the above results 
slightly modify. The main physical reason is that the neutrino plasma receives a small 
contribution from the e + e~ annihilations, and its final energy density is a little bit higher 
than in the standard case. Neutrinos with higher momenta are more heated since, in the 
range of energies we are interested in, weak interactions get stronger with rising energy. 
This causes a momentum-dependent distortion in the neutrino spectra from the pure 
equilibrium Fermi-Dirac shape. 

There have been a number of papers which considered the effects of non-instantaneous 
neutrino decoupling [||, f|, ||, || [7|, |[ |9j . The first papers used integrated quantities to esti- 
mate that the neutrino energy density increases by a factor 1% || |], [5|, while subsequent 
works H |7| made some momentum-dependent calculations assuming Boltzmann statistics 
for neutrinos and other approximations. The full numerical computation of the evolution 
of the neutrino distribution functions without approximations requires the numerical so- 
lution of the Boltzmann equations, as done in [BL The most accurate results for the 
evolution of the neutrino distortion until complete decoupling are given in Ref. O. For 
what concerns the energy stored in the relic neutrinos, this study gives a temperature 
ratio after decoupling T 7 /7j, = 1.3991 , which indeed shows that neutrinos also share a 
small part of the energy transfer from the e + e~ annihilations. The method adopted in 



this computation has however been questioned in Ref. [[HJ (see also JTT|] for a reply to 
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this work). One of the aims of our work is to verify the results of Ref. ||. We have 
also performed a full numerical calculation, although employing a different method. Our 
analysis confirms the results of Ref. || . 

The main goal of our study is to extend the previous analyses to the case in which 
a large asymmetry between neutrinos and antineutrinos is present. If this asymmetry is 
of order one, one says that neutrinos are Fermi-degenerate. Big Bang Nucleosynthesis 
(BBN) forces the baryonic asymmetry 7]^ = {n^ — 77,5) /n 7 to be very tiny (of the order 
10~ 10 — 10 -9 ) [0, and the observed electric neutrality of the Universe translates these 
bounds on electrons too. However, the existence of a large asymmetry (even of order one) 
in the neutrino sector is an open possibility which has drawn much attention in the past. 

Actually, there exist some realistic theoretical models that can generate very large 
values of the lepton asymmetry 7] v . From a particle physics point of view, a lepton 
asymmetry can be generated by an Affleck-Dine mechanism [13| without producing a 



large baryon asymmetry (see Refs. ]14|, |H| [16], [17|, 0), or even by active-sterile neutrino 



oscillations JHJ after the electro-weak phase transition. In general, the asymmetry is 
expected to be different for each neutrino family. 



Important bounds on r\ v come again from BBN (see [21] for a review of non-standard 
BBN). The asymmetry on electronic neutrinos is the most constrained, since it has a 
direct influence on the weak interactions between protons and neutrons and thus affects 
the final 4 He abundance. Asymmetry on the neutrinos of the two other families is instead 
less bounded, since its only effect is to increase the expansion rate of the universe. The 
simultaneous presence of an asymmetry in the three neutrino families can still lead to a 
successful prediction of the abundances of light elements (fully degenerate BBN), provided 
that their values are suitably chosen. This last possibility usually requires the presence of 
a baryon asymmetry somewhat larger than the one allowed by standard nucleosynthesis 
(that is with non-degenerate neutrinos). Indeed, a high neutrino asymmetry can be 
considered the simplest option to save BBN predictions in case a high baryon asymmetry 



is required. The very recent Boomerang p2j and Maxima p3f results on the acoustic peaks 



of the Cosmic Microwave Background Radiation (CMBR) seem to favour ||24| , P5[, g6[ P7fl 
high values for rj b , which may be in conflict with the bounds coming from standard BBN. 
In particular, in Ref. |27| it is argued that the presence of a neutrino asymmetry may 
contribute to both (i) improve the fits of the Boomerang data and (ii) render the high rjb 
needed compatible with (degenerate) BBN. 

Finally, a relic neutrino asymmetry, delaying the matter domination stage, can have 
very significant effects on the matter power spectrum. In particular, it can suppress the 
power at small with respect to large scales, thus making the predictions of Cold Dark 
Matter models compatible with observations. 

Motivated by these considerations, in this work we extend the study of non-equilibrium 
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effects on the relic neutrino spectra to the case of non- vanishing neutrino chemical poten- 
tials (degenerate neutrinos). The paper is organized as follows. In section 2 we review the 
effect of neutrino chemical potentials on BBN and other cosmological observables. The 
set of equations ruling the evolution of neutrino distributions and the method adopted 
to get the solution are presented in section 3. Numerical details are given in section 4, 
whereas section 5 contains our results. We discuss in section 6 some phenomenological 
implications of the non-equilibrium effects. Finally, in section 7 we give our conclusions. 



2 Cosmological implications of degenerate neutrinos 

The existence of a relic neutrino asymmetry would have important cosmological conse- 
quences mainly on Primordial Nucleosynthesis and the CMBR and matter power spectra. 

Degenerate neutrinos influence BBN in two distinct ways. The first one is connected 
to the increase of the energy density of the primordial plasma due to a non-vanishing 
neutrino asymmetry. One can introduce an effective number of relativistic neutrinos N v 
which, for degenerate neutrinos, depends on the chemical potential of any neutrino flavour. 
The quantity N v is defined, in fact, from the total neutrino (plus antineutrino) energy 
density, p„, through the relation 

which, for three massless neutrinos with chemical potentials \i Va and in the equilibrium 
case, becomes 



N v = 3+ 



OL = e,fJ,,T 



(2) 



where £ a = \i Va jT v are the degeneracy parameters. Note that N v does not depend on the 
sign of £ a . For what concerns BBN, N v > 3 leads to a higher neutron to proton ratio 
since it favours an earlier decoupling of nucleons. This produces an increase in the final 
production of both 4 He (because there are more neutrons available when the nuclei form), 
D and 3 He (since there is less time to destroy them). 

A second important effect on BBN predictions is induced by the asymmetry for v e 
only, since they directly participate in the n <-> p weak processes. An excess of electronic 
neutrinos over antineutrinos (£ e > 0), in fact, enhances the n — > p conversion rate with 
respect to the inverse process, so reducing the neutron to proton ratio. In addition, the 
initial n/p ratio diminishes by a factor exp(— £ e ), and this further reduces the number of 
neutrons available at the onset of BBN. As a consequence, the observed abundances of 
light elements can be achieved also in this case, but with a value of Q b h 2 significantly 
higher than for £ e = 0. 



4 



When more than one neutrino species is degenerate, both the above effects combine 
and, as a result, one can observe particular combinations of values of £ a for which the 
predictions of degenerate BBN are still in good agreement with the observational data on 
the abundances of primordial elements. 

An exhaustive analysis of degenerate BBN was performed in Ref. |28| (see refs. therein 



for previous works), while more recently some aspects have been studied in [29, 30, 31 
Neglecting the non-electron neutrino chemical potentials, the authors of [^HJ find the limits 
—0.06 < £ e < 0.14. Instead, for fully degenerate BBN and requiring that the Universe had 
a sufficiently extended period of matter domination, the neutrino degeneracy parameters 
lie in the wider ranges |28| 



- 0.06 < £ e < 1.1 and { \^\ % ™ ** ^ ^ ° ' f = ^ vice Versa ) (3) 

[ IC^rl < 5.6 for £ M = £ T ^ 

since, as mentioned above, at least for 4 He, the effect of a positive £ e can be compensated 
by the contribution to N u coming from £ M)T . 



In a very recent work [32], the BBN bounds on the neutrino degeneracies have been 
re-analyzed. The input parameters of degenerate BBN are £ e , N u (£ a ) and r] W = lO 10 ?]^. 
Using an updated code, which includes all relevant physical effects that influence the 4 He 
abundance up to 0.1%, one can perform a likelihood analysis of compatibility between 
theoretical predictions and experimental data. This analysis yields contour levels which 
are surfaces in the three-dimensional space of these parameters. However, one should take 
into account other non BBN constraints for reducing the ranges of the parameters. For 
example, considering the upper bound on the radiation density present at recombination 



coming from CMBR data, ref. |33| obtains a 2a limit N u < 13. As far as £ e and r/io 
are concerned, the analysis is limited to the range — 1 4- 1 and 1 4 30, respectively. The 
maximum for these functions is found for |53] 

£ e = 0.06 , iV, = 3.43 , »7io = 5 , (4) 

for a low value of the D abundance, and 

£ e = 0.35 , iV, = 13 , r/io=4.20 , (5) 



in the high D case (see figures 19 and 20 of [32 for the 95% exclusion plots for the £ e and 



N u parameters for different values of 7710) . For low D the allowed range for 7710 is 3.3 -=-9.9, 
while for high D we have 1.1 -j- 5.8. 

Primordial nucleosynthesis is not the only framework which can be substantially 
affected by a non-zero relic neutrino degeneracy. For example, in the recent past, several 
papers have considered the imprint of this asymmetry on the power spectra of CMBR 
anisotropies and matter density. It has been found |34], |55|, that the asymmetry 
boosts the amplitude of the first CMBR peak, shifts the peaks to larger multipoles, and 
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suppresses small scale matter fluctuations (see |J7|] for a previous work). All these effects 
are the consequences of increasing the neutrino energy density, which delays the epoch of 
matter-radiation equality, and from their analysis one can put constraints on the neutrino 
degeneracy |34], |36|. Among the most recent observations, some interesting features have 



emerged from the new Boomerang results [22], which show a quite puzzling suppression 
of the second acoustic peak with respect to the first one. A possible explanation to this 
problem is provided by the simultaneous presence of both high baryonic and high neutrino 



asymmetries, as shown in Ref. [2?|. The possible consequences of the neutrino degeneracy 
on the future CMBR measurements (Planck) are discussed in Ref. . 



As far as the matter power spectrum is concerned, it is shown in Ref. |36| that 
a non-vanishing neutrino degeneracy can suppress the small scale fluctuations, so as to 
agree with the observations even in a pure Cold Dark Matter scenario. It is also shown 



in Ref. |36| that the suppression at small scales is particularly efficient if the degenerate 
neutrinos are massive, because free-streaming of non-relativistic neutrinos is enhanced 
when their average momentum is boosted by the chemical potential. 

The existence of a relic lepton asymmetry enhances the contribution of massive 
neutrinos to the present energy content of the Universe |B8|, Actually it has 

been shown that even the smallest neutrino mass suggested by Super-Kamiokande data 



on atmospheric neutrinos |4(| could be extracted from CMBR anisotropy and large-scale 
structure data by the future Planck satellite and Sloan Digital Sky Survey, provided that 



a large neutrino asymmetry exists 41 



Finally, it is worthwhile observing that there could be other implications of the 
presence of degenerate relic neutrinos, that include the explanation of the ultra-high 
energy cosmic rays, beyond the Greisen-Zatsepin-Kuzmin cut-off of about 5 x 10 19 eV. 
These cosmic rays would be produced by the protons from the annihilation of ultra-high 
energy neutrinos on the relic degenerate neutrinos in the galactic halo (more abundant 
than in the standard case) at energies close to the Z-resonance. 



3 The dynamics of neutrino distributions 

At the time of neutrino decoupling, the evolution of their distributions functions is de- 
scribed by a set of Boltzmann equations where the collisional terms are due to the weak 
interactions of neutrinos with the primordial plasma. Since the baryonic component is 
much smaller than the leptonic one, we can safely neglect its presence. In this case, the 
whole set of relevant reactions are those reported in Ref. ||, where a complete study 
of the evolution of neutrino distributions for the non-degenerate case, namely vanishing 
chemical potentials, is performed. In this section we describe a method which allows us 
to extend the results of Ref. || to the degenerate case. 
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Following the notations of Ref. we choose as time variable x = m e R, and 
comoving momentum y = p u R, where R is the universe scale factor and m e is the electron 
mass. We also define the rescaled photon temperature aS Z = TryR. 

At sufficiently high temperatures neutrinos can be considered in thermal equilibrium 
with the 7, e plasma through weak interactions (we always consider temperatures below 
the muon mass^). Thus they are described by Fermi-Dirac distributions with the same 
temperature of the electromagnetic plasma. However, at lower temperature one expects 
a different evolution for the distributions of electronic neutrinos and antineutrinos, which 
also experience charged current interactions due to the presence of in the thermal bath, 
with respect to the neutrinos of the other families. Thus, in the following, we will assume 
identical distributions for muon and tau neutrinos (antineutrinos), 

fu„ = fur = fv x { X ,V) > /*„ = fih = fu x ( X ,V) ) ( 6 ) 

so restricting the unknown neutrino distributions to f Ue , fp e , f Ux and f Bx only. 

In the temperature range we are interested in, electrons and positrons are kept in 
thermodynamical equilibrium with photons by fast electromagnetic interactions. Thus, 
they are distributed according to the Fermi function | 

f e -(x,y,z) = f e +(x,y,z) = 1 — — . (7) 

exp(\/y z + x z /z) + 1 

In order to get the time evolution of neutrino distributions, f Ua (x, y), and the rescaled 
photon temperature, z(x), one must solve the following set of differential equations 

d , „ 1 / 

—p(x) = -(p-3P 
ax x v 

^f^{ x ,y) = I y a {fu e Ju e Ju x Ju x ] , with V a = U e , v e , u x , v x . (9) 

In equation (§), which states the conservation of the total energy density, p and P are 
the dimensionless energy density and pressure of the primordial plasma, respectively, 

\ 4 / \ 4 



P = P[ — ) . p = p — , (10) 

\m e J \m e J 

and the index m in the r.h.s. reminds that only massive components in the plasma con- 
tribute. From (H), by using the expressions 

P7 = Y / ' 



Un the range of ^ v we arc interested in the vv annihilation gets out of equilibrium at temperatures 
well below the muon mass ||^, Q . 

*We neglect the completely irrelevant asymmetry since it is expected to be £ e ~ r\b < 10 -9 . 



7 



3 it 2 Jo y/x 2 + y 2 exp(v / y 2 + x 2 / z) + 1 ' 
1 r°° 

P u = o 2 / ^ ^ L/W^y) +Mx,y) + 2f Ux (x,y) + 2f Px (x,y)} , (14) 



2tt 2 Jo 

we get the equation for the evolution of z{x 



<iz z 4z 3 jo y c?x dx dx dx 

dx x 2 2tt a 

-F^x z) + F 2 (x z) + — 

z z 15 

where the functions are given by 



(15) 



OJ 2 + T 2> 



I a 2 exp(y^ . , , 
Fl(T) = («xp(^ + T») + l)» • (16) 



cxp 



(Vo; 2 + t 2 ) 



F 2 (t) = duu 4 . (17) 

V ' Jo {exp{VuJ 2 + t 2 ) + l) 2 V ; 

From eq. fll5l) , neglecting the terms proportional to the derivative of neutrino distributions, 
one gets the asymptotic value z® q = (11/4) 1 / 3 which represents the ratio between the 
photon and neutrino temperatures after the complete annihilation of e + e~ pairs. In the 
presence of neutrino chemical potentials, the behaviour of z(x) is in general different. 
However, the final value of z(x) is always lower than (11/4) 1 / 3 , showing that also the 
neutrino plasma is slightly heated by the e + e~ annihilations (mathematically, this can be 
seen from the fact that the non-equilibrium contributions involving neutrino distributions 
in the r.h.s of ( |T5*D are of negative definite sign). 

In the set of Boltzmann equations @, I Va represents the collisional integral for the 
single neutrino species v a , and is a functional of all neutrino and e ± distributions. At the 
time of neutrino decoupling, the plasma density was low enough that, in the expression 
of I Va , one can safely consider only two-body weak reactions 1 + 2^3 + 4 with v a = 1, 



Ii/ a [fu e j fu e j fv x i fv x \ — 2 jjj J 



d 3 p2 d 3 p3 d 3 p4 



\3 



2 * reasons 7 2 E 2 (2^ 2 E 3 (2^ 2E 4 (2n) 
x (2vr) 4 ^(p 1 +p 2 -p3-P4) F[/l,/2,/ 3 ,/4]|M 12 ^34| 2 (18) 

where F = f 3 f 4 (1 - fi) (1 - f 2 ) - h /a (1 - fa) (1 - h) is the statistical factor, and 
^12^34 is the process amplitude. In Ref. the complete list of relevant processes and 
corresponding squared amplitudes are reported, and it is shown that using the 5-function 
some of the integrals can be analytically performed, reducing I Voi to a two-dimensional 
integral. 
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In order to solve eqs. (Bj) and QI5D, instead of using a discretization in momentum 



space (discrete values for y, see Refs. || ||), we employ a method based on orthonormal 
polynomials^. According to this general technique, we rewrite the four unknown neutrino 
distribution functions as 

U a {x,y)= 1 (l + Sf Va (x,y)) , (19) 

e y U + 1 

where £ a is the neutrino degeneracy parameter. The initial conditions on f Ua {x,y) are 
fixed by observing that it is possible to find a starting value for the evolution parameter 
x = Xi n , such that for temperatures larger than the corresponding T in , the neutrino 
distributions can be safely assumed to be the equilibrium ones. This is envisaged by the 
Fermi factor in the r.h.s of fljjD, provided a vanishing 8f„ a for x = X{ n . In this case the 
initial conditions depend on the two independent input parameters £ e and £ x only, being 
the antineutrino equilibrium distributions characterized by the corresponding opposite 
chemical potentials. We take x in = m e /(10 MeV) as the starting value for the evolution 
parameter (as in 0). Solving the equilibrium part of eq. (Kf), the initial z is found to be 
z in = z(x in ) = 1.00006. 



According to (|T9|) , the function 5f Ua parameterizes the departure from equilibrium 
and can be expanded in terms of a set of polynomials, P"(y), as 

oo 

*AJ*,v) = £<*?(*) J?(v) > (2°) 

i=0 

where Pf{y) are constructed, with the standard Gram-Schmidt orthonormalization pro- 
cedure, and the requirement that they are orthonormal with respect to the Fermi function 
weight, 

f idrrT^M^)^ • < 21 > 

Note that each set of polynomials depends on the neutrino degeneracy parameter £ a 
through (PT|). By substituting ( ffiD in eqs. (BJ), we can rewrite them as 

d_ 

dx 

with u a = u e , u e , u x , v x and % — 0, 1, .... From the above considerations, the initial condi- 



1 r°° 

a?(x) = — j o d Vl P?( yi ) I Va [f, e ,f, e ,f, x ,h x ] , (22) 



tions for (22) are af(xi n ) = 0. 



4 Numerical details 



In order to solve eqs. (|T5|) and (E^), we have to truncate the infinite series in eq. (EU|) at 
a term % = m. The choice of m is driven by the requested accuracy, which we take to 

§Note that an expansion of the non-equilibrium distortions in momenta was also discussed in ref. |E3 
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be of the order 1% in the neutrino distortions. One can make an estimate of the error in 
approximating eq. ( ^Oj) by comparing the results for two subsequent values m and m + 1. 
We have verified that the requested accuracy can be obtained retaining the coefficients 
until m — 3, 

W,»)^EfliW^) • (23) 

8=0 

If the expression of each polynomial is considered, 

P l a (y) = E 6g 2/ j , (24) 

one can rewrite Eq. ( p3[) as 

8f„ a (x,y)~j2<$(z)y i , (25) 



i=0 

where the coefficients c" are given by 

c?(x) = E h% a«(x) . (26) 

3=0 

The evolution equations for the n = 4(m + 1) + 1 unknown functions af and z(x) are 
solved with an integrator for stiff equations. This is implemented in the Fortran code by 
calling a NAG routine for stiff equations, which uses Backward Differentiation Formulas 
with Newton's method and an adaptive step-size. In order to speed up the evaluation of 
the r.h.s. of Eq. (|i"5|) , we made a fit of the functions F\ and F 2 with a precision better than 
1%. On the other side, as far as the 4(m + 1) equations (^) are concerned, we needed to 
compute the tridimensional integral on the r.h.s. of eq. (0). At this aim we used a routine 
implementing the Korobov-Conroy number theoretical method, repeating the integration 
twice for calculating the errors, and made a check on the results by comparing them with 
the ones obtained with other integrators. 

Actually, one of the equations, that is the energy conservation law flT5|), is not stiff. 
Thus in order to check our results for some values of the neutrino chemical potentials, we 
used two versions of the code: in the first one we solved straightforwardly eqs. (fl5|) and 
( p2f) as previously described, whereas in the second one we first separately considered the 
equilibrium component in z(x), which satisfies 

x 

d ~Fi{x/z) 

—F 1 (xz)+F 2 (xz) + — 
z z 15 

In Figure[I|, the quantity z eq is plotted versus x. It is independent on the neutrino chemical 
potentials and asymptotically gives the value z® q = (11/4) 1 / 3 . Then, the remaining non- 
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equilibrium correction Az( ) satisfies the differential equation 

1 r A 3 ( d U , df De df Ux df Dx \ 



d . Az 3 Jo \ dx dx dx dx I , . 

* *V \ 777, 2 ** ' ( ' 



which has a vanishing r.h.s. once the neutrinos completely decouple and their distribution 
functions get frozen (df /dx = 0). The two methods have produced very well coinciding 
results, and this is a check of the good performance of the NAG integrators. 



5 Results 

In order to compare our results with the analysis of Ref . || , we have first solved Eqs. fll5|) 
and (p2[) for vanishing neutrino chemical potentials (£ e = £ x = 0). In Figures and ^| the 



distortion coefficients c\ and c?, defined in eq. (pq) , as a function of the time variable x are 
reported. Note that the antineutrino distributions coincide with the neutrino ones since 
we have vanishing chemical potentials. From the plots one can see that the coefficient 
evolution ends at x ~ 10, which means T 7 = m e z(x)/x ~ 0.07 MeV. After this value the 
coefficients cf reach their asymptotical values and the neutrino distributions can be well 
defined as thermodynamically decoupled from the electromagnetic plasma. Note that in 
both cases, the value of the coefficient is always much smaller than the others, justifying 
the truncation of the expansion in (p5) (in reasonable agreement with the estimate of 0). 



In Figure £| the total distortions 5f Ue and Sf Vx as a function of x are reported, 
for three values of neutrino momentum y. Note that such distortions increase with the 
neutrino energy and, as expected, they are more relevant for electron neutrinos, due to 
the fact that only these interact with e^ through charge currents. This is clear from 
Figure where we can see that the final distortion Sf Ua is larger for v e than for v x and 
is an increasing function of neutrino energy. The contribution of the distortion to the 
differential energy density is given in Figure || where one can see that the maximum is for 
y ~ 5. All the results obtained here for the non-degenerate case are in perfect agreement 
with the ones reported in 0. 

In Figure [7] we plot the evolution of Az/z eq in % with x as obtained from the 
solution of eq. (ESI) for different choices of neutrino degeneracy parameters. In all cases, it 
is a negative decreasing function approaching a constant value when neutrinos decouple. 
Thus, from this plot, one can see that the neutrinos are completely decoupled (i.e. the 
distortions 5f are frozen) at x ~ 3. One also notices that Az (the energy transfer from 
e + e _ to neutrinos) slightly diminishes when the neutrino degeneracies increase. This is a 

^Note that the bump at x ~ 3 — 10 is only due to the different values of z and z eq at the late stages 
of their evolution, and would not appear if z(x) and z eq (x) were plotted separately. 
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consequence of the fact that the rate of the e + e annihilation channel into neutrinos is 
reduced in the presence of neutrino degeneracy. 

In Figures ||| and |] we instead report the evolution of the total distortions for all 
neutrino and antineutrino species as a function of x for given values of the neutrino 
degeneracy parameters (£ e = 0.5, £ x = 1) and for three values of neutrino momentum 
y. Comparing these results with the ones for the non-degenerate case, we observe that 
the distortions are enhanced for antineutrinos and depressed for neutrinos. The different 
behaviour between v and v is also evident if one considers an effective Fermi distribution 
with a degeneracy parameter £ e jy as a function of the momentum. It is easy to check 
that, at first order in the neutrino distortion, £ e ff is given by 

C ff (x, y) ~ T + (1 + exp(r - y))8fu(x, y) 
e eff (x, y) ~ -C + (1 + exp(-r - y))6fp(x, y) (29) 



In Figures |TU] and |TT] we plot the absolute values of £ e ff for a fixed momentum y = 5. 
The total final distortions and their effect on the differential energy density are plotted 
in Figures ^ and One can see from the last plot that the maxima of the neutrino 
distortions are displaced from y ~ 5 by an approximate factor ±£„ a . 

We have solved the differential equations flip] ) and ( ]2"2| ) for < £j, e < 0.5 and 
< £ Vx < 1 obtaining, for all these values, the final distribution functions f„ a (y)- The 
general expression for this quantity as a function of the neutrino degeneracy parameters, 
in the above ranges, has been fitted with the following form: 



% M = y-L I i + U, k) a ^ V 1 exp {B? & + C? & ) . ( 30 ) 

' i" 1 \ i=0 j,k=0 



The coefficients of these fits are reported in Tables We have checked that in order 
to obtain the extension of fl30f) for negative values of £ Q is a good approximation to 
exchange the neutrino distribution functions with those of antineutrinos. An example is 



given in Figures [14] and [15] for the cases £ x = ±0.6, both calculated solving the differential 



equations (|15D and (|22| ) 



6 Phenomenological implications of the neutrino dis- 
tortions 

There are two cosmological scenarios where the non-equilibrium effects on the neutrino 
spectra could in principle manifest: Primordial Nucleosynthesis and the spectrum of 
CMBR anisotropics. 
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We have calculated the change in the production of primordial elements when the 
evolution of the neutrino distortions is included. The non-equilibrium neutrino distribu- 
tions produce three different effects. The first one is to increase the total energy density of 
the Universe during BBN and this produces a larger 4 He abundance. The inclusion of 8f Va 
also modifies the weak n <-> p processes, essentially increasing the rates, which has as an 
overall effect the destruction of some neutrons (and a smaller 4 He abundance), since they 
are more abundant than protons. Finally the evolution law of the photon temperature is 
also changed to take into account the small energy transfer from e + e~ to neutrinos. Since 
the final value of the photon temperature is slightly smaller than in the standard case, 
this last effect leads also to a less effective production of 4 He. In order to facilitate the 
comparison of our results with those of previous works, we have incorporated the distor- 



tions 5f Ua (x,y) to the classic BBN code of Kawano f44 |. We have found that the three 
effects of the neutrino distortions almost cancel, and the final change in the mass fraction 
(Y) of 4 He is only at the 10~ 4 level, in good agreement with previous analyses || |9|, |45| , 



which is undetectable given the current observational uncertainties. In Figure we show 
AF as a function of the neutrino degeneracies (the changes in the other elements are even 
smaller) . 

The non-equilibrium effects that we have calculated will slightly modify the relic 
value of the neutrino energy density. We show in Figure [T7| the increase of the neutrino 



energy density, parameterized in terms of an excess in the effective number of neutrinos 
N„ as in eq. (0), 

This small increase in p v can have a significant effect on the CMBR. The epoch of matter- 
radiation equality is very sensitive to the radiation energy density of the Universe, and 
the slight heating of neutrinos changes the predicted spectrum of CMBR anisotropics. 
Even if the effect is too small to be noticed with the present data, it has to be taken into 
account when determining the other cosmological parameters [ 46[ . It has been found that 



the small correction to the present neutrino energy density is marginally detectable if the 
anisotropy and polarization of the CMBR is measured with the expected precision for the 
satellite mission Planck ESI . 



7 Conclusions 

We have calculated the exact kinetic evolution of neutrinos in the early Universe un- 
til complete decoupling, extending previous analyses to the case in which a large neu- 
trino asymmetry exists. We numerically solved the Boltzmann equations for the neutrino 
distribution functions, finding the momentum-dependent corrections to the equilibrium 
spectrum of neutrinos. At the same time, we found the final value of the ratio between 
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photon and neutrino temperatures, which is slightly smaller than in the standard case due 
to small energy transfer from the electromagnetic plasma to neutrinos due to e + e _ — > vv 
processes. 

Our results are in nice agreement with those of refs. || [Tl| for the non-degenerate 
case, obtained with a different method of solving numerically the kinetic equations. In- 
stead of a discretization of the distribution function in momenta, we employed an expan- 
sion in momenta of the non-equilibrium part, truncated at order y 3 for an accuracy of 1%. 
Let us present our results for this case: z = 1.39905 for the final photon temperature, 
8p Ve /Pu e = 0.953 % and 8p v / p v = 0.399 % for the corrections to the energy density of 
electron and muon/tau neutrinos, respectively. 

When the neutrino degeneracies are non-zero, we found that the distortions in the 
neutrino and antineutrino momentum spectra are different, and larger for antineutrinos 
(neutrinos) for positive (negative) asymmetries. We have obtained the final distribution 
functions in the range of neutrino degeneracies < £„ e < 0.5 and < £ Vx < 1, and 
presented the corresponding fits in Tables [L]-^. 

We have finally discussed two cosmological scenarios where the existence of distor- 
tions to the spectra of cosmic neutrinos could have phenomenological implications. For 
what concerns Primordial Nucleosynthesis, we have found that the various effects of the 
distortions almost cancel, leading to increases in the final mass fraction of 4 He of the order 
10~ 4 with a small dependence on the neutrino degeneracies, well below the present obser- 
vational uncertainties, in agreement with previous studies. Finally the distortions make 
a small positive contribution Sp u to the relic value of the neutrino energy density, which 
slightly decreases in the presence of neutrino degeneracies. This contribution increases 
the radiation content of the Universe and can have an effect on the CMBR that must be 
taken into account when determining the other cosmological parameters [}fE|| . 
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Table 1: Values of the coefficients 10 2 -A" e (j, k). 
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Table 2: Values of the coefficients 10 2 -Af e (j, k). 
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Table 3: Values of the coefficients 10 2 -A^(j, k). 
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Table 4: Values of the coefficients 10 2 -A^(j, k). 
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Table 5: Values of the coefficients B" a and C" a . 
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Figure 1: Evolution of the quantity z eq (x), defined in eq. (|27|). The asymptotical value 
represents the well-known value = (11/4) 
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Figure 2: Evolution of the coefficients function of x (non-degenerate case). 
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Figure 4: The evolution of the total distortion 8f Va (25), as a function of x (non-degenerate 
case) for three values of neutrino momentum y. From bottom to top: y = 3, y = 5 and 
y = 7. 
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Figure 6: Final distortion of the differential energy density, as a function of momentum 
(non-degenerate case). 
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Figure 7: The evolution of Az/z eq from eqs. (|27|) and (|28|), as a function of x for four 
choices of the neutrino degeneracy parameters £ e , 




Figure 8: The evolution of the total distortion Sf Va , as a function of x (£ e = 0.5, £ x = 1) 
for three values of neutrino momentum y. From bottom to top: y = 3, y = 5 and y = 7. 
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Figure 9: Same as previous figure, for the antineutrino distortions. 



0.51 




J 0.495 - 
0.49 - 
0.485 - 



0.48 1 ■ ■ ■ 1 

0.1 1 10 

x 

Figure 10: Evolution of the absolute value of the effective degeneracy parameter £ (if 
the distribution function is written with an equilibrium form) for electron neutrinos and 
antineutrinos with momentum y = 5 (£ e = 0.5, ^ = 1). 
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Figure 11: Same as previous figure, for muon or tau neutrinos and antineutrinos. 
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Figure 16: Effect of the neutrino distortion over the primordial abundance of 4 He as a 
function of the neutrino degeneracies £ e and £ x . The contours indicate equal values of AY 
in units 10~ 4 . 
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Figure 17: Effect of the neutrino distortion over the relic neutrino energy density, para- 
meterized in terms of AN U , as a function of the neutrino degeneracies £ a . 
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